Higher Charmonia 
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This paper gives results for the spectrum, all allowed El radiative partial widths (and some 
important Ml widths) and all open-charm strong decay amplitudes of all 40 cc states expected 
up to the mass of the 4S multiplet, just above 4.4 GeV. The spectrum and radiative widths are 
evaluated using two models, the relativized Godfrey-Isgur model and a nonrelativistic potential 
model. The electromagnetic transitions are evaluated using Coulomb plus linear plus smeared 
hyperfine wavefunctions, both in a nonrelativistic potential model and in the Godfrey-Isgur model. 
The open-flavor strong decay amplitudes are determined assuming harmonic oscillator wavefunctions 
and the 3 Po decay model. This work is intended to motivate future experimental studies of higher- 
mass charmonia, and may be useful for the analysis of high-statistics data sets to be accumulated 
by the BES, CLEO and GSI facilities. 
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I. INTRODUCTION 

Since its discovery in 1974 0,0, the charmonium sys- 
tem has become the prototypical 'hydrogen atom' of me- 
son spectroscopy 0, 0, IS 19- The experimentally clear 
spectrum of relatively narrow states below the open- 
charm DD threshold of 3.73 GeV can be identified with 
the IS, IP and 2S cc levels predicted by potential mod- 
els, which incorporate a color Coulomb term at short 
distances and a linear scalar confining term at large dis- 
tances. Spin-dependent interquark forces are evident in 
the splittings of states within these multiplcts, and the 
observed splittings are consistent with the predictions of 
a one gluon exchange (OGE) Breit-Fermi Hamiltonian, 
combined with a linear scalar confining interaction. Dis- 
cussions of the theoretical importance and experimental 
status of heavy quarkonium, including recent experimen- 
tal results for charmonium, have been given by Quigg 
Galik the CERN quarkonium working group @ , Seth 
EEllll and Swarnicki 0. 

Recently there has been a resurgence of interest in 
charmonium, due to the realization that B factories can 
contribute to the study of the missin g cc states , and 
to high-statistics experiments at BES |l5j and CLEO 0] 
and the planned GSI pp facility 0J ■ 
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The possibility of contributions from B factories was 
dramatically illustrated by the recent discovery of the 
long missing 2 1 So r]' c state by the Belle Collaboration [Isj . 
which has since been confirmed by BABAR Il9| . and has 
also been observed by CLEO in 77 collisions |2J| . 

Additional interest in cc spectroscopy has followed the 
discovery of the remarkable X(3872) by Belle [2l| and 
CDF [23 in B decays to J / 'ybTr + n~ ; assuming that this 
is a real resonance rather than a threshold effect, the 
X(3872) is presumably either a DD* charmed meson 
molecule [23|, |2J, |25j or a narrow J = 2 D-wave cc state 
[US- Very recent observations of the X(3872) in 7 J ftp 
and wJ I ip by Belle support a 1 ++ DD* molecule assign- 
ment [28ll29| . 

There has also been experimental activity in the spin- 
singlet P-wave sector, with recent reports of the obser- 
vation of the elusive l^Pi h c state by CLEO 0, l3fj| . 
Finally, the surprisingly large cross sections for dou- 
ble charmonium production in e + e _ reported by Belle 
[31TI32I l33l l34l | suggest that it may be possible to study 
C = (+) cc states in e + e~ without using the higher-order 
0{a A ) two-photon annihilation process. 

One open topic of great current interest in cc spec- 
troscopy is the search for the ^(l 3 !^) and ?y C 2(l 1 E>2) 
states, which are expected to be quite narrow due to the 
absence of open-charm decay modes. 

A second topic is the Lorentz nature of confinement; in 
pure cc models this is tested by the multiplet splittings of 
orbitally-excited cc states. For example, with pure scalar 
confinement as is normally assumed there is no spin-spin 
hyperfine interaction at 0(v 2 /c 2 ), so the masses of spin- 
singlets (such as the x Pi h c ) are degenerate with the 
corresponding triplet c.o.g. (center of gravity); here this 
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is the 3 P,| c.o.g., at 3525 MeV. In the original Cornell 
model [35J it was assumed that confinement acts as the 
time component of a Lorentz vector, which lifts the de- 
generacy of the h c and the 3 Pj c.o.g. Another possibil- 
ity is that confinement may be a more complicated mix 
of scalar and timelike vector j^| . Of course these sim- 
ple potential model considerations may be complicated 
by mass shifts due to other effects, such as couplings to 
open-flavor channels |27| . 



A third topic is the search for exotica such as hybrids; 
the level of mixing between conventional quarkonium and 
hybrid basis states falls rapidly with increasing quark 
mass, which suggests that nonexotic hybrids may be more 
easily distinguished from conventional quarkonia in char- 
monium than in the light quark sectors. Since LGT (lat- 
tice gauge theory) predicts that the lightest cc hybrids lie 
near 4.4 GeV [371 HE EM EH , there is a strong incentive 
to establish the "background" spectrum of conventional 
cc states up to and somewhat beyond this mass. 

A final topic of current interest is the importance of 
mixing between quark model qq basis states and two- 
meson continua, which has been cited as a possible reason 
for the low masses of the recently discovered D s j states 
[4ll |42| . The effects of "unquenching the quark model" 
by including meson loops can presumably be studied ef- 
fectively in the cc system, in which the experimental spec- 
trum of states is relatively unambiguous. The success of 
the qq quark model is surprising, in view of the probable 
importance of corrections to the valence approximation; 
the range of validity of the naive "quenched" qq quark 
model is an interesting and open question 0] . 

Motivated by this revived interest in cc spectroscopy, 
we have carried out a theoretical study of the expected 
properties of charmonium states, notably the poorly un- 
derstood higher-mass cc levels above DD threshold. Two 
variants of potential models are used in this study, a con- 
ventional nonrclativistic model based on the Schrodinger 
equation with a Coulomb plus linear potential, and the 
Godfrey-Isgur relativized potential model. We give re- 
sults for all states in the multiplets 1-4S, 1-3P, 1-2D, 
1-2F and 1G, comprising 40 cc resonances in total. Pre- 
dictions are given for quantities which are likely to be of 
greatest experimental interest, which are the spectrum 
of states, El (and some Ml) electromagnetic transition 
rates, and strong partial and total widths of states above 
open-charm threshold. 

Similar results for many of the electromagnetic tran- 
sition rates have recently been reported by Ebert et al. 
p4|. The £ + £~ leptonic and two-photon widths arc not 
discussed in detail here, as they have been considered ex- 
tensively elsewhere; see for example ^5j and 0, 0, Esf 
and references cited therein. 



II. SPECTRUM 
A. Nonrelativistic potential model 

As a minimal model of the charmonium system we 
use a nonrelativistic potential model, with wavefunc- 
tions determined by the Schrodinger equation with a con- 
ventional quarkonium potential. We use the standard 
color Coulomb plus linear scalar form, and also include 
a Gaussian-smeared contact hyperfine interaction in the 
zeroth-order potential. The central potential is 
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where 8 a (r) = (a/^/n) 3 er a r . The four parameters 
(a s , b, m c , cr) are determined by fitting the spectrum. 

The spin-spin contact hyperfine interaction is one of 
the spin-dependent terms predicted by one gluon ex- 
change (OGE) forces. The contact form oc 5(x) is actu- 
ally an artifact of an 0(u 2 /c 2 ) expansion of the T-matrix 



so replacing it by an interaction with a range 1/a 
comparable to 1 /m c is not an unwarranted modification. 

We treat the remaining spin-dependent terms as mass 
shifts using leading-order perturbation theory. These are 
the OGE spin-orbit and tensor and a longer-ranged in- 
verted spin-orbit term, which arises from the assumed 
Lorentz scalar confinement. These are explicitly 
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The spin-orbit operator is diagonal in a |J,L, S) basis, 
with the matrix elements (L • S) = [J(J + 1) — (L(L + 
1) — S(S + l)]/2. The tensor operator T has nonvanishing 
diagonal matrix elements only between L > spin-triplet 
states, which are 
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For experimental input we use the masses of the 11 
reasonably well established cc states, which are given 
in Table [I] (rounded to 1 MeV). The parameters that 
follow from fitting these masses are (a s , 6, m c , a) = 
(0.5461,0.1425 GeV 2 , 1.4794 GeV, 1.0946 GeV). Given 
these values, we can predict the masses and matrix el- 
ements of the currently unknown cc states; Tableland 
Fig- n show the predicted spectrum. 



B. Godfrey-Isgur relativized potential model 

The Godfrey-Isgur model is a "relativized" extension 
of the nonrelativistic model of the previous section. This 
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model assumes a relativistic dispersion relation for the 
quark kinetic energy, a QCD-motivated running coupling 
a s (r), a flavor-dependent potential smearing parameter 
a, and replaces factors of quark mass with quark kinetic 
energy. Details of the model and the method of solu- 
tion may be found in Ref . |Hflj . The Hamiltonian consists 
of a relativistic kinetic term and a generalized quark- 
antiquark potential, 

H = H +V q q(p,r) (4) 

where 

Ho = y/p% + m\ + v /5f + m§ . (5) 

Just as in the nonrelativistic model, the quark- 
antiquark potential V q q(p, r) assumed here incorporates 
the Lorentz vector one gluon exchange interaction at 
short distances and a Lorentz scalar linear confining in- 
teraction. To first order in (v q /c) 2 , Y qq -(p,r) reduces 
to the standard non- relativistic result given by Eqns.(l) 
and (2) (with a s replaced by a running coupling con- 
stant, aJr)). The full set of model parameters is given 
in Ref. 50] . Note that the string tension and quark mass 
(b = 0.18 GeV 2 and m c = 1.628 GeV) are significantly 
larger than the values used in our nonrelativistic model. 

One important aspect of this model is that it gives 
reasonably accurate results for the spectrum and ma- 
trix elements of quarkonia of all u, d, s, c, b quark flavors, 
whereas the nonrelativistic model of the previous section 
is only fitted to the cc system. 

C. Discussion 



a constraint on the spectrum of higher-L cc states gen- 
erally. This is unfortunately a difficult problem, since 
these states are not easily accessible experimentally. As 
we shall see, one possibility is to produce the 3 F2 X2 cc 
state, which is formed in El radiative transitions from the 
-0(4159) and decays dominantly to DD and DD*. (For 
simplicity in this paper we abbreviate the final state DD 
as "DD", the state DD* +DD* as "DD*", and so forth.) 
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FIG. 1: Predicted and observed spectrum of charmonium 
states (TablcQJ. The solid lines are experiment, and the bro- 
ken lines are theory (NR model left, GI right). Spin triplet 
levels are dashed, and spin singlets are dotted. The DD open- 
charm threshold at 3.73 GeV is also shown. 



The spectra predicted by the NR and GI models (Ta- 
ble [I] and Fig. ^) are quite similar for S- and P-wave 
states, largely because of the constraints provided by the 
experimental cc candidates for these multiplets. We note 
in passing that these potential model results are very 
similar to the most recent predictions of the charmonium 
spectrum from LGT pM loll l52j ] . At higher L we have 
only the L=2 l 3 Di and 2 3 Di states 0(3770) and ^(4159) 
to constrain the models, and the predicted mean D-wave 
multiplet masses differ by ca. 50 MeV. For L > 2 the ab- 
sence of experimental states allows a relatively large scat- 
ter of predicted mean masses, which differ by as much as 
« 100 MeV in the 1G multiplet. (The splittings within 
higher L multiplets in contrast are rather similar.) The 
mean multiplet masses predicted by the two models dif- 
fer largely because of the values assumed for the string 
tension b, which is 0.18 GeV 2 in the GI model but is a 
rather smaller 0.1425 GeV 2 in the NR model. Identifica- 
tion of any L=3 or L=4 cc state would be very useful as 



III. RADIATIVE TRANSITIONS 
A. El transitions 

Radiative transitions of higher-mass charmonium 
states are of interest largely because they provide one of 
the few pathways between cc states with different quan- 
tum numbers. Since typical El radiative partial widths 
of charmonia are 10s to 100s of keV, corresponding to 
significant branching fractions of ~ 10 -3 to 10~ 2 , large 
event samples of radially excited S-wave states produced 
in e + e~ annihilation could be used to identify radially 
excited P-wave states, which are not otherwise easily pro- 
duced. Similarly, the El radiative transition of the nom- 
inally 2 3 Di -0(4159) can be used to produce an F-wave 
cc state; this multiplet would likely be very difficult to 
reach using other mechanisms. 

We evaluate these El radiative partial widths using 
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r E1 (n *+% - n' 2S '+ 1 L' J , + 7) = \ C fi S ss , e\ a | | r \ fr) | 2 E 7 — 7-jy , (6) 



where e c = 2/3 is the c-quark charge in units of |e|, a is 
the fine-structure constant, E 7 is the final photon energy, 
e[ cc ' ) is the total energy of the final cc state, Mj 00 ' is the 
mass of the initial cc state, the spatial matrix element 
{if)f\ r I ipi) involves the initial and final radial wavefunc- 
tions, and the angular matrix element Cfi is 

C/^max^, L')(2J' + l)|j^| 2 . (7) 

This is the result of Ref . [ 53j , except for our inclusion of 

the relativistic phase space factor of E^ - 1 , which is 

usually not far from unity. (The GI results do not include 
this phase space factor.) We evaluate these radiative par- 
tial widths in both the NR and GI models. For the NR 
model the matrix elements (n' 2S +1 L'j,| r |n 2S+1 Lj) were 
evaluated using the Coulomb plus linear plus smeared hy- 
perfine wavefunctions of the potential model described in 
Sec.IIA, and for the GI model they were evaluated using 
the wavefunctions of Ref. 50] . Since the masses predicted 
for unknown states differ in the two models, the assumed 
photon energy E 7 differs as well; these photon energies 
are given in the El and Ml transition tables (Tables ITTl 
I1X|) together with the radiative partial widths. 

Some El transitions that are of special importance for 



(See the previous El formula for definitions.) The GI 
Ml radiative rates do not incorporate the E^ /M^ 00 ' 
phase space factor, but do include a jo(fcr/2) recoil factor. 
We quote NR results both with and without this recoil 
factor. The photon energies E 7 depend on the model in 
most cases, since we have assumed theoretical masses for 
unknown states. 

As with the El rates, we quote Ml radiative widths in 
both the NR potential model of Sec.IIA and the GI model 
[50| of Sec.IIB. Although the Ml decay rates involve a 
relatively simple off-diagonal matrix element of the mag- 
netic moment operator, we note that our predicted rates 
show considerable variation with the model assumptions, 
and agreement with the two known Ml rates is not good. 
In the off-diagonal (2S — ► IS) transition ip' — > jr] c this 
is because the zeroth-order Ml matrix element vanishes 



I 

the study of higher charmonium states are discussed in 
the text. Transitions from initial 1" cc states are of 
greatest interest in this regard, since these can be stud- 
ied with high statistics at e + e~ machines. These can 
provide access to the spin-triplet members of the 2P and 
3P multiplets in particular, starting from the -0(4040) 
and -0(4415). El radiative transitions may also be useful 
in identifying the narrow ID 3 _ and 2~ cc states, since 
they are all predicted to have large partial widths (ca. 
300 keV) to the IP xj an d h c states. 



B. Ml transitions 

Although Ml rates are typically rather weaker than 
El rates, they are nonetheless interesting because they 
may allow access to spin-singlet states that are very dif- 
ficult to produce otherwise. It is also interesting that 
the known Ml rates show serious disagreement between 
theory and experiment. This is in part due to the fact 
that Ml transitions between different spatial multiplets, 
such as tp 1 — » 7ry c (2S — > IS), are nonzero only due to 
small relativistic corrections to a vanishing lowest-order 
Ml matrix element. 

The Ml radiative partial widths are evaluated using 



(8) 

I 

due to orthogonality of the spatial wavefunctions, and 
the nonzero predicted rate results from rather model- 
dependent corrections to the wavefunctions and recoil 
effects. Evidently the corrections we have included do 
not accurately predict the observed partial width; for 
this reason the other predicted Ml rates between differ- 
ent spatial multiplets are suspect, and in any case are 
evidently strongly dependent on recoil factors. Better 
experimental data will be very important for improving 
our description of these apparently simple but evidently 
poorly understood Ml radiative transitions. 



I 

T M1 (n 2S ^L 3 -> n' 2S '+ 1 L; i , + 7 ) = I 2 ^±± 6 W * s , s , ±1 e 2 ± | A) ? ^ 
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IV. OPEN-FLAVOR STRONG DECAYS 
A. The Decay Model 

The dominant strong decays (when allowed by phase 
space) are transitions to open-flavor final states. In 
these open-flavor decays the initial cc meson decays 
through production of a light qq quark- antiquark pair 
(q = u, d, s), followed by separation into two open-charm 
mesons. Remarkably, the QCD mechanism underlying 
this dominant decay process is still poorly understood. In 
quark model calculations this decay process is modelled 
by a simple phenomenological qq pair production ampli- 
tude. The qq pair is usually assumed to be produced 
with vacuum (0 ++ ) quantum numbers, and variants of 
the decay model make different assumptions regarding 
the spatial dependence of the pair production amplitude 
relative to the initial cc pair. The simplest of these mod- 
els is the 3 P model, originally introduced by Micu [54| . 
which assumes that the new qq pair is produced with 
vacuum quantum numbers ( 3 Po) by a spatially constant 
pair-production amplitude 7. 

LeYaouanc et al. subsequently applied the 3 P model 
to meson and baryon |5fil open flavor strong de- 
cays in a series of publications in the 1970s. They also 
evaluated strong decay partial widths of the three cc 
states_j/>(3770), 0(4040) and ^(4415) in the 3 P model 
|58l l59j : the relation between this early work and our 
contribution will be discussed. 

This decay model, which has since been applied ex- 
tensively to the decays of light mesons and baryons, was 
originally adopted largely because of the success in pre- 
dicting the D/S amplitude ratio in the decay b\ — ► ujt: 
[55l l&OL loll l63| . Another stringent test of strong decay 
models is provided by the very tight limit on the spin 
singlet to spin singlets transition 7^(1670) — > bin from 
the VES Collaboration IHHi, 

5^(1670)^6^ < 1-9 • 10- 3 , 97.7% c.l. (9) 

This branching fraction is predicted to be zero in the 
3 Po model, but would not necessarily be negligible in 
a different decay model or if final state interactions are 
important. (Final state interactions combined with the 
3 Po model allow two-stage transitions such as tt2 — *• pit — > 
&i7r, although the direct decay 712 —> biir is forbidden in 
the 3 P model.) 

Recent variants of the 3 Po model modify the pair pro- 
duction vertex 65] or modulate the spatial dependence of 
the pair production amplitude to simulate a gluonic flux 
tube |60| . (The latter is the "flux-tube decay model" , 
which gives very similar predictions to the 3 Po model in 
practice.) 

Another class of decay models assumes that the pair 
production amplitude transforms as the time component 
of a Lorentz vector. The Cornell group used a decay 
model of this type for charmonium [35j . This model ap- 
pears to describe the partial and total widths of some 
charmonium states reasonably well, although it is known 



to disagree with experimental amplitude ratios in the 
light quark sector |61j. 

In this work we employ a formalism that is equiva- 
lent to the original constant-amplitude version of the 3 Po 
decay model, although we have simplified the calcula- 
tions by deriving momentum-space Feynman rules |6l| 
instead of using the real-space convolution integrals of 
LeYaouanc et al. |55j . 

In our formalism the 3 Po model describes decay matrix 
elements using a qq pair production Hamiltonian which 
is the nonrelativistic limit of 

Hdecay = 7 ^Z^ m q / ^ X g ? . (10) 

q J 

Here ^ q is a Dirac quark field, m q is the constituent quark 
mass, and 7 is the dimensionless qq pair production am- 
plitude, which is fitted to data. The decay arises from 
the matrix element of this model Hamiltonian between 
an initial meson state |A) and a final meson pair |BC), 
which is nonzero due to the b'cP term. Here we will use 
SHO (simple harmonic oscillator) wavefunctions for the 
mesons, with a universal width parameter /?. The pair- 
production strength parameter 7 is found to be roughly 
flavor-independent in light meson decays involving pair 
production of uu, dd and ss pairs. In our recent exten- 
sive studies of light (u, d, s) meson decays |(01|6(||67j we 
found that the value 

7 = 0.4 (11) 

gives a reasonably accurate description of the overall 
scale of decay widths. 

B. Charmonium Strong Decays 

1. Previous studies 

The open-flavor decay amplitudes of the three cc states 
■0(3770), 0(4040) and ^(4415) were evaluated in the 3 P 
model by LeYaouanc et al. in the late 1970s [58l l59j . 
Their calculations are generally quite similar to the ap- 
proach used here, although we note three important dif- 
ferences: 

1) LeYaouanc et al. fitted the pair production strength 
7 and wavefunction length scale R to the charmonium 
decay data, whereas we used the same pair production 
amplitude 7 as in light meson decays, and our SHO wave- 
function length scale (3~ l is taken from cc potential mod- 
els. 

2) We use the correct reduced mass coordinates for the 
unsymmetric open-charm meson wavefunctions, whereas 
LeYaouanc et al. assumed symmetric SHO wavefunctions 
for all mesons. This approximation proves to be quite 
important numerically, since m c 3> m Sl m u ^. 

3) We consider all energetically allowed open-flavor de- 
cay modes. LeYaouanc et al. did not consider some S+P 
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decays of the -0(4415), nor did they consider the impor- 
tant J p = 1 + D meson singlet-triplet mixing angle. (The 
masses of the P-wave charmed mesons were unknown at 
that time.) 

Although we have improved on the earlier calculations 
of LeYaouanc et al, and evaluate a much more exten- 
sive set of decay amplitudes, we do concur with their 
important conclusion that the DD mode of the 0(4040) 
is naturally suppressed given a 3 3 Si cc assignment, and 
also that S+P modes are important for the ip(4A15). 

Other more recent applications of the 3 Po decay model 
to charmonium include a study of the strong decays of 
various cc X(3872) candidates to DD [2^|, and an exam- 
ination of the hypothesis that the ■0(4040) and 0(4160) 
are linear superpositions of a cc hybrid and a conven- 
tional 0(3S) cc state [H. 

The only other detailed studies of the open-flavor 
strong decays of charmonia of which we are aware are 
the early work by the Cornell group |35| , and recent more 
detaile d ap plications of this model to other charmonium 
states 27] . The decay model of the Cornell group as- 
sumes a pair creation operator which is the time compo- 
nent of a Lorentz vector (j°), rather than the Lorentz 
scalar assumed by the 3 Po model. As we noted pre- 
viously, this j model does not agree well with some 
light meson decay amplitude ratios. In the original Cor- 
nell model decay calculations only the combined chan- 
nel contributions to R were derived, and only channels 
containing pairs of S-wave open-charm mesons were in- 
cluded. The very interesting relative amplitudes in the 
D*D* channel have not been evaluated in this model. 



2. This study 

Tables IXIXVI present the partial widths and strong 
decay amplitudes we find for all kinematically allowed 
open-flavor decay modes of all the charmonium states 
listed in Table |U We have evaluated these in the 3 Po 
model, assuming SHO wavefunctions with a width pa- 
rameter of /3 = 0.5 GeV and a pair production ampli- 
tude of 7 = 0.4. SHO wavefunctions have the advantage 
that the decay amplitudes can be determined analyti- 
cally, and we have found that the numerical results are 
usually not strongly dependent on the details of the spa- 
tial wavefunctions [g£l l6ll 16^. Fof , unless they are near 
a node. This width parameter was chosen by comparing 
the overlap of NR and GI potential model wavefunctions 
with SHO wavefunctions, which was found to be largest 
for both charmonium and open-charm mesons for a (3 
near 0.5 GeV. 

As a test of the accuracy of these parameters, espe- 
cially the value assumed for the qq pair production ampli- 
tude 7 (which is taken from light meson decays), in Fig. [21 
we compare the predicted total widths for (3 — 0.5 GeV 
and variable 7 to the values observed for the four known 
cc states above open-charm threshold, -0(3770), 0(4040), 
0(4159) and 0(4415). Although there is some scatter 



in the value of 7 specified by these widths, a choice of 
7 = 0.35 is evidently near optimum, and yields an av- 
erage error of only 29%, which is very reasonable for 
this phcnomcnological model. In the decay tables we 
will quote numerical results for our standard light-meson 
value of 7 = 0.4, which is rather close to this value. A 
subsequent change in 7 will simply scale the widths as 
7 2 - 
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FIG. 2: Values of the 3 Po decay model pair production 
strength 7 implied by the experimental total widths of the 
known higher-mass vector states -0(3770), 0(4040), 0(4159) 
and 0(4415). This calculation assumes pure cc spectroscopic 
states, respectively l 3 Di, 3 3 Si, 2 3 Di and 4 3 Si. The boxes 
show current PDG experimental width uncertainties. 



For this paper the strong decay amplitudes were first 
evaluated analytically, following which numerical values 
were determined given our parameters (3 = 0.5 GeV, 7 = 
0.4, and our constituent quark masses m Ui d = 0.33 GeV, 
m s = 0.55 GeV and m c = 1.5 GeV. (Quark mass ratios 
are required to specify the final open-charm meson wave- 
functions.) We determine the decay kinematics using ex- 
ternal meson masses taken from the Particle Data Group 
|64| (charge avera ged where appropriate) and from re- 
cent Belle results [71|. These masses are Md = 1.867 
GeV, M D * = 2.008 GeV, M D * = 2.308 GeV [Hi, 
M Dl (narrow)= 2.425 GeV, M D ?(broad)= 2.427 GeV 
M D ; = 2.459 GeV, M Da = 1.968 GeV, M D . = 2.112 
GeV, M D - o = 2.317 GeV, M Dsl = 2.459 GeV. The FO- 
CUS collaboration estimate a somewhat higher mass 
for the scalar Dg, however their result is complicated by 
the presence of both scalar and broad axial vector (D^) 
contributions [73| . 

The J p = 1 + axial vector cn and cs mesons Di and 
D' x are assumed to be coherent superpositions of quark 
model spin-singlet and spin-triplet states, 



|Di) = +cos(0)| 1 P 1 ) 
|Di) =-^(0)1^0 



-sin^)! 3 ?!) 

■COs(0)| 3 Pi) 



(12) 



We define this singlet-triplet mixing angle in the LS cou- 
pling scheme, in accord with the conventions of Ref . 67] . 



This reference also discusses other conventions for this 
mixing angle that have appeared in the literature. In the 
heavy-quark limit we expect to find a "magic" mixing 
angle, due to the quark mass dependence of the spin- 
orbit and tensor terms, which is 6 — 35.3°(— 54.7°) if the 
expectation of the heavy-quark spin-orbit interaction is 
negative (positive) [74|. In the following we assume the 
first case, which is supported by the reported widths of 
the 1 + P-wave cn mesons. We note however that finite 
quark mass effects and mixing induced by higher-order 
Fock states can substantially modify this mixing angle, 
so it should more generally be treated as a free parame- 
ter. We will discuss the dependence of our strong decay 
amplitudes on this mixing angle for a few sensitive cases. 

V. DISCUSSION OF cc STATES 
A. Known states above 3.73 GeV 




-90 -45 45 90 



6 [deg.] 

FIG. 3: The total width of the 0(3770) as a function of the 
3 Di - 2 3 Si mixing angle 6. The experimental la limits are 
shown as horizontal dashed lines. 



The four known cc states above DD threshold, -0(3770) , 
0(4040), -0(4159) and 0(4415) are of special interest be- 
cause they are easily produced at e + e~ machines. Ac- 
cordingly we will first discuss our predicted strong decay 
amplitudes and widths for these states. This will be fol- 
lowed by a more general discussion of cc strong decays 
by multiplet. 



1. 0(3770) 

The 0(3770) is generally assumed to be the 3 D i c c 
state, perhaps with a significant 2 3 Si component |75ll7q . 
(This additional component can explain the leptonic 
width, which is much larger than expected for a pure 
3 Di cc state 0.) 

For a pure 3 Di state at the 0(3770) mass we predict a 
DD width of 43 MeV with our parameters, which is rather 
larger than the experimental value of 23.6±2.7 MeV. The 
partial width of a mixed 2S-D state 



10(3770)) = +cos(6»)| 3 Di) +sin(6»)|2 3 Si 



(13) 



is shown in Fig. [31 as a function of the mixing angle 8. 
In this simple model, fitting the experimental 0(3770) 
width requires a mixing angle of 9 = —17.4° ± 2.5°. 
Assuming that the leptonic widths scale as the S-wave 
component squared, this mixing angle predicts an e + e _ 
partial width ratio of 



r e+e - (0(3770)) 



(0(3686)) 



= 0.10 ±0.03 



(14) 



thy. 



It is interesting that this is consistent with experiment, 
(0(3770)) 



r e + e - 



r e + e 



(0(3686)) 



= 0.12 ±0.02 



(15) 



expt. 



so the small 0(3770) total width may indeed be an effect 
of 3 Di - 2 3 Si mixing. 



2. 0(4040) 

Next in mass is the 0(4040), which is a very inter- 
esting case for the study of strong decays. Four open- 
charm modes are energetically allowed, DD, DD*, D*D* 
and D S D S , although D*D* has little phase space. There 
is some experimental evidence for the three nonstrange 
modes from Mark I at SLAC [77j . Remarkably, the re- 
ported relative branching fractions (scaled by p 3 ) show 
a very strong preference for D* final states, D*°D*° >> 
D°D*° >> D°D°. This motivated suggestions that the 
0(4040) might be a D*D* molecule [73, US The 
Mark I results were 



B 

P 



3 ( D*°D*° 



D°D*° : D°D° ) 
128 ±40 : 4.0 ±0.8 : 0.2 ±0.1 



(16) 



This is much less striking once the p 3 factors are restored, 
which gives 



B( D*°D*° : D°D*° : D°D° ) = 
1±0.31 : 0.95 ±0.19 : 0.12 ±0.06 



(17) 



where we have used the D*°D*° mode to set the scale. 
These relative branching fractions should be compared 
to our theoretical 0(4040) branching fractions, not the 
more often quoted B/p 3 ratios. 

Comparison with Tablelxlshows that the predicted 3 Po 
model partial and total widths for a 3 3 Si 0(4040) state 
are actually in accord with the Mark I results. The total 
width is predicted to be 74 MeV, which is reasonably 
close to the PDG experimental average of 52 ± 10 MeV. 
Seth [8l| estimates a rather larger width of 88 ± 5 MeV 
for this state from Crystal Ball and new BES data, which 
is closer to our predicted width. 

The DD* and D*D* branching fractions are predicted 
to be approximately equal, consistent with Mark I. 
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The DD mode is predicted to be much smaller, despite 
its larger phase space, as it is accidentally near a node 
in the 3 Po model decay amplitude. (This was previously 
noted by LeYaouanc et al. |5fi|.) This set of predictions 
constitutes a very nontrivial success of the 3 Po decay 
model. The actual width to the suppressed DD mode is 
probably a few MeV; the 0.1 MeV quoted in the table 
depends strongly on the location of the node, and vary- 
ing the SHO width parameter f3 by ±10% suggests that 
a DD partial width of a few MeV is plausible. This is 
again consistent with the Mark I result. 

The D*D* mode is especially interesting for strong 
decay studies, since there are three independent decay 
amplitudes for 1~ cc -> D*D*, 1 P 1 , 5 Pi and 5 Fi. If 
the 0(4040) is a pure S-wave cc state, we predict a zero 
5 Fi D*D* decay amplitude, and the ratio of the nonzero 
amplitudes is 5 Pi/ 1 Pi = —2^5 (independent of the ra- 
dial wavefunction) . As we shall see when we discuss the 
'0(4159), very different D*D* amplitude ratios are pre- 
dicted for a 3 Di initial cc state. 

The unobserved mode D S D S is predicted to have a 
branching fraction of about 11%; this prediction is fortu- 
nately not especially sensitive to the location of the decay 
amplitude node for 0(4040) to two pseudoscalars. This 
branching fraction is of special interest because it deter- 
mines the event rates available for studies of D s weak 
decays at e + e~ colliders such as CLEO. 

Finally, we note that a large sample of -0(4040) events 
might be used to access the 2P cc multiplet, since the 
El radiative branching fractions for -0(4040) — > 7Xj are 
expected to be ca. 10~ 3 (Table |B}. 



3. 0(4159) 

Comparison of the mass of the 1~~ -0(4159) with po- 
tential model predictions (Fig. ^| immediately suggests 
a 2 3 Di cc assignment. There may also be a significant S- 
wave cc component, since the -0(4159) has a much larger 
e + e~ width than one would expect for a pure D-wave cc 
state 0- (This was also the case for the 0(3770).) 

There are five open-flavor decay modes available to 
a 1 cc vector at this mass, DD, DD*, D*D*, D S D S 
and D S D*. The predicted decay amplitudes and partial 
widths for a 2 3 Di 0(4159) are given in Table IXlTl The 
theoretical total width of 74 MeV is in good agreement 
with the experimental value of 78 ± 20 MeV. Seth |8lj 
quotes a similar width of 107 ± 8 MeV for this state from 
Crystal Ball and recent BES data. 

The leading mode is predicted to be D*D* with a 
branching fraction of « 50%, followed by comparable 
DD and D S D* modes (both « 20%), and a somewhat 
weaker D S D S (« 10%). The DD* branching fraction is 
predicted to be very small, since it is suppressed by a 
decay amplitude node near the physical point. 

The mode D*D* is again especially interesting, due to 
the three decay amplitudes allowed for this final state. 
(The other modes have only single amplitudes.) For a 



pure D-wave cc assignment the ratio of the two D*D* P- 
wave amplitudes is independent of the radial wavefunc- 
tion, and is 5 Pi/ 1 Pi = -l/y/E. The 5 Fi amplitude is 
predicted to be the largest given this 2 3 Di assignment, 
whereas it is zero for an S-wave cc state. Clearly, a de- 
termination of these D*D* decay amplitude ratios would 
be an extremely interesting test of decay models. Exper- 
imentally, to date nothing has been reported regarding 
the exclusive hadronic decay modes of the -0(4159). 



4. 0(4415) 

The final cc resonance known above DD threshold is 
the 1 -0(4415). Again, the mass of this resonance 
relative to potential model predictions suggests a cc as- 
signment, in this case 4 3 Si. Of course this identification 
requires independent confirmation, since cc hybrids are 
predicted to first appear near this mass by LGT simula- 
tions [U IH, |Zj| , and the lightest hybrid multiplet in 
the flux tube model includes a 1"" state |82^ . 

Ten open-charm strong decay modes are allowed for 
the 0(4415), seven with cn meson final states (n — u, d), 
and three with cs. The predicted branching fractions of a 
4 3 Si cc state at this mass are quite characteristic, and (if 
the 3 Po decay model is accurate) may be useful in con- 
firming this assignment. As a note of caution, a 4S state 
has three radial nodes, and some of the smaller predicted 
branching fractions are sensitive to the locations of the 
nodes. 

The predicted total width of a 4 3 Si cc meson at 
this mass with our parameters is 77 MeV, somewhat 
larger than the experimental PDG average width of 
43 ± 15 MeV. Seth 81] notes however that the width of 
this state in Crystal Ball and recent BES data is rather 
larger, and quotes an average of 119 ± 15 MeV, on the 
opposite side of the 3 Po decay model prediction. 

The largest exclusive mode is predicted to be the S+P 
combination DDi, where Di is the narrower of the two 
1 + cn axial mesons near 2.42-2.43 GeV. Since the Di is 
rather narrow (r « 20-30 MeV) and decays dominantly 
to D*7r, there should be a strong 0(4415) signal in DD*7r 
final states. Although both S-wave and D-wave DDi final 
states are allowed, in the 3 Po model the HQET (heavy 
quark effective theory) Di mixing angle 6 is just the value 
needed to give a zero S-wave 0(4415) — + DDi amplitude. 
Thus we have the striking prediction that the dominant 
0(4415) decay mode is DDi, in D-wave rather than S- 
wave. 

The second-largest 0(4415) branching fraction is pre- 
dicted to be another S+P mode, DD?,. The D?, is also 
moderately narrow, so this isobar can also be isolated 
from the observed final state. (The PDG quotes a neu- 
tral D£ total width of T = 23 ± 5 MeV, but Belle 
fn\ and FOCUS [8^ find somewhat broader values of 
r = 45.6 ± 8.0 MeV and T = 38.7 ± 5.3 ± 2.9 MeV re- 
spectively; see |U for a recent discussion. The D?, has 
significant branching fractions to both D*n and Dtt, so 
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the DD*, mode of the 0(4415) should be observable in 
both DDtt and DD*tt. 

A final important -0(4415) mode is predicted to be 
D*D*, which should be comparable in strength to DD*,. 
As noted in the previous discussions of 0(4040) and 
-0(4159) decays, D*D* is an especially interesting decay 
mode because it has three amplitudes, and the 5 Pi/ 1 Pi 
amplitude ratio is independent of the radial wavefunction 
for pure S-wave or D-wave cc states. If the 0(4415) is in- 
deed an S-wave (4 3 Si) cc state to a good approximation, 
we expect this ratio to be 5 Pi/ 1 Pi = — 2\/5, and the 5 Fi 
amplitude should be zero. 

It is interesting to note that 0(4415) decays may pro- 
vide access to the recently discovered D* (2317). Al- 
though the channel D*D* (2317) has a threshold of 
4429 MeV, 14 MeV above the nominal mass of the 
0(4415), the width of the 0(4415) and the fact that 
the decay 0(4415) -^D*D* (2317) is purely S-wave im- 
plies that one may observe significant D* (2317) pro- 
duction just above threshold, near E cm sa 4435 MeV. 
This is illustrated in Fig. in which we show theoreti- 
cal 3 Po decay model partial widths to D*D* (2317) and 
D s D s i(2459) as functions of the 4 3 Si cc mass. (This cal- 
culation assumes pure csD* (2317) and D s i(2459) states, 
and should be modified accordingly if they have signifi- 
cant non-cs components.) Unfortunately, 0(4415) decays 
are not expected to be similarly effective in producing the 
D s i(2459), because the assumed "magic-mixed" HQET 
state |D B i(2459)) = y/T/3 | 3 Pi(cs)) + | x Pi(cs)) 
predicts a vanishing 0(4415) — > D s D s i S-wave decay am- 
plitude. 




4.42 4.43 4.44 4.45 4.46 4.47 4.48 4.49 
M(\|/(4S)) [GeV] 

FIG. 4: Partial widths predicted for 0(4 3 Si) -> D*D* (2317) 
and D s D s i(2459) as a function of the assumed 0(4 3 Si) mass. 



B. 3S and 4S states 

The two unknown states in the 3S and 4S multiplets 
are the 3 1 So and 4 1 So pseudoscalars. To evaluate their 
strong decays we have assigned potential model masses 
to these states, specifically 4043 MeV and 4384 MeV. 
The resulting total widths are predicted to be moderate, 
80 MeV and 61 MeV respectively, and both states should 
be observable in DD* and D*D* (see Table 0. No other 
open-charm strong decay modes are allowed for a 3 1 Sq 
7y c (4043). Although these modes are also important for 
the 4 1 So ?7c(4384), its largest branching fraction is pre- 
dicted to be to the S+P combination DD*, (w 40%). 

These states may be observable in 77 and hadronic 
production (through gg fusion), analogous to the rj c and 
rf c . Ml transitions from the higher vectors to these states 
are unfortunately predicted to be rather weak, since they 
are either hindered or have little phase space. (See Ta- 
ble El) 



C. 2P states 

The mean 2P multiplet mass is predicted to be near 
3.95 GeV. There is an interesting disagreement between 
the NR and GI models regarding the scale of mass split- 
tings within this multiplet, as well as in the 3P states. 
(See Table H| and Fig. [J) 

There are few open-flavor strong modes available to the 
2P states (Table 1X1}) . One result of this restricted phase 
space is the prediction of a fairly narrow 2 3 Po xo(3852) 
(assuming the NR model mass) that decays only to DD, 
with a total width of just 30 MeV. This small width is 
due in part to a decay amplitude node, and a variation of 
the total width by a factor of two could easily be accom- 
modated. The two axial states 2 3 Pi xi(3925) and 2^1 
/i c (3934) can only decay to DD*, and both have nonzero 
S- and D-wave decay amplitudes; their D/S ratios are 
—0.18 and +0.46 respectively. Observation of the D/S 
ratio in a DD* enhancement could be a useful check of 
the assumption of a resonant 1 + cc contribution. 

The relatively narrow 2 3 Po state xo(3852) can be pro- 
duced in 77, as can its somewhat wider 2 3 P2 X2(3972) 
partner. Additional possible production mechanisms in- 
clude gg fusion (for the 2 3 Po and 2 3 P2), El radiative 
transitions from higher 1 — - cc states (to all 2P xj states), 
and B decays. The El branching fractions of the 0(4040) 
and 0(4159) to 2P xj states are predicted to be suffi- 
ciently large (ca. 10~ 3 , see Tables ITIlV|) to allow identifi- 
cation of these states at BES and CLEO. 

Although no 2P cc state has been clearly established 
experimentally, there are two recent reports of enhance- 
ments which may be due to states in this multiplet. The 
Belle Collaboration has reported evidence of an enhance- 
ment in ujj/tp with a mass and width of 3941 ±11 MeV 
and 92 ± 24 MeV [8J|; this is compatible with expec- 
tations for 2P xj cc states. (Although 2P cc states 
should strongly favor open-charm strong decay modes, 
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inelastic FSIs (final state interactions) will allow weaker 
transitions to closed-charm final states such as coJ/ip. 
These FSIs should be most important in S-wave, as in 
the virtual second-order processes Xo,2(3940) — > D*D* — > 
lo J ftp.) A search for evidence of this state in open-charm 
final states, with much larger branching fractions than 
to uJ/ip, will be a crucial test of a charmonium assign- 
ment. Assuming that the enhancement is indeed due to 
a 2P xj state, if it is the ++ state it will populate only 
DD, if 1 ++ only DD*, and if 2 ++ comparable branching 
fractions to DD and DD* are expected. 

There is a second report of a possible 2P state, also 
from the Belle Collaboration [33|, in double charmo- 
nium production. An enhancement is seen in e + e~ — > 
J/V> DD* at a mass of 3940 ± 12 MeV, with a limit on 
the total width of < 96 MeV (90% c.L). As this enhance- 
ment is not seen in tvj/ip, it has been suggested that 
it is distinct from the 3940 MeV signal discussed above. 
Consideration of possible 2P candidates (Table 1X1(1 sug- 
gests that this enhancement may be due to the Xci(2P), 
which is predicted to have a similar mass of 3925 MeV. 
(The predicted total width however is a somewhat larger 
165 MeV.) 



D. 3P states 

The 3P states have an expected mean multiplet mass of 
about 4.3 GeV. As with the 2P states there is a significant 
difference between the NR and GI models in the scale of 
mass splittings predicted within the multiplet. This leads 
for example to a 90 MeV difference between the models 
for the predicted mass of the 3 3 Po scalar. 

Many open-charm strong decay channels are open at 
the 3P mass scale. The individual decay amplitudes 
however are typically somewhat smaller than for the 2P 
states, due in part to nodes in the decay amplitudes. This 
has the interesting result that the predicted mean total 
width of 3P states is smaller than for 2P states (58 MeV 
for 3P versus 91 MeV for 2P). 

The predicted branching fraction to the mode D*D* 
is large for all four 3P states, and is the largest for all 
except the spin-singlet 3 : Pi /i c (4279); in this case the 
leading mode is DD*>. For all 3P states except the 3 3 Pi 
the D*D* decay mode has multiple amplitudes, which 
are predicted to be comparable in strength. An exper- 
imental determination of any of these amplitude ratios 
would provide an interesting test of decay models. The 
decays of the tensor 3 3 P 2 X2(4317) to DDi and DD^ are 
also interesting, in that the prediction that the broad D^ 
mode dominates is a test of the axial-D mixing angle. 

As with the 2P states, the most important experimen- 
tal problem is to find an effective production mechanism. 
Again, 77 and gg fusion can produce the 2 ++ and ++ 
states (here 3 3 P2 and 3 3 Po), although there should be 
moderate suppression of the production amplitudes rela- 
tive to the IP ground states due to the higher masses and 
smaller short-distance cc wavef unctions. El transitions 



from the ?/>(4415) could lead to the identification of the 
3P states, since radiative partial widths of 50-100 keV 
are expected to the 3 3 P2 and 3 3 Pi; these correspond to 
branching fractions of ca. 10~ 3 . Production of the 3 3 Po 
state however is expected to be suppressed, due to a ra- 
diative decay amplitude node. 



E. ID states 

Since the 3 D2 and 1 D2 states do not have allowed 
open-flavor decay modes, and we have already consid- 
ered 3 Di decays in the section on the -0(3770), we re- 
strict our discussion to the 3 D3. This currently unknown 
state is predicted to have a mass of 3806 MeV in the 
NR potential model. It is expected to be quite narrow 
simply because it has little phase space to the only open- 
flavor mode DD, which has an F-wave centrifugal barrier; 
with our parameters the predicted total strong width is 
just 0.5 MeV. The radiative width for the El transition 
-0 3 (38O6) -> 7x2 is predicted to be about 0.3 MeV (Ta- 
blelVjl. so the total width of this 3 D3 state should be near 
1 MeV, with comparable strong and radiative branching 
fractions. 

Finding an effective production mechanism for this in- 
teresting state is again a crucial problem. As it has 
C = (— ), it cannot be made in 77 or gg fusion. One possi- 
bility is that the 3 D3 may be observable in El transitions 
from the radial 2 3 P2 tensor, since the radiative partial 
width is predicted to be a relatively large « 100 keV. 
Similarly, the narrow 3 D2 and 1 D2 states may be observ- 
able in El transitions from 2P states, assuming that a 
production mechanism can be found for 2P states. All 
three of these states can also be produced in B-decays 
0. 



F. 2D states 

The states in the 2D multiplet are predicted to be es- 
sentially degenerate in both the NR and GI potential 
models. The 2 3 Di candidate -0(4159) suggests a 2D mul- 
tiplet mass near 4.16 GeV, which is consistent with po- 
tential model expectations. Since the ^(4159) was dis- 
cussed previously, here we will consider only the unknown 
states 2 3 D 3 , 2 3 D 2 and 2 1 D 2 . 

These states are predicted to be rather broad, with the- 
oretical total widths ranging from 78 MeV for the 2 3 Di 
(consistent with the -0(4159)) to 148 MeV for the 2 3 D 3 . 
The energetically allowed open-flavor modes are all of 
S+S type, since the first S+P mode is DDq, at a nominal 
mass of 4175 MeV. The decays of 2D states are domi- 
nated by DD, DD* and D*D*, with smaller contributions 
from the charm-strange final states D S D S and D S D*. 

If these states (besides the VK4159)) can be produced, 
there are several interesting tests of their strong decay 
amplitudes. The 2 3 D 3 ?^ 3 (4167), for example, is pre- 
dicted to have D*D* as the dominant mode, with both 
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1 F 3 and 5 F 3 amplitudes present, in the ratio 5 F 3 / 1 F 3 = 
~^2AjZ- The remaining 2D states have many multi- 
amplitude decays, in which the higher-L partial waves 
have large or dominant amplitudes. (Compare the P- 
and F- waves in the decays in Table IXTTl ) 



G. IF states 

The four states in the IF multiplet are predicted to 
be almost degenerate in both the NR and GI models, 
although the models disagree somewhat regarding this 
mass. In the NR potential model the states are expected 
near 4025 MeV, whereas in the GI model the expected 
mass is near 4095 MeV. This discrepancy is mainly due 
to the different values assumed for the string tension. 

It is unfortunate that production amplitudes for these 
higher-L states are expected to be quite weak, since the 
IF multiplet is predicted to contain a rather narrow state, 
the 3 F 4 X4 ( Table IXiTT|) . Assuming our NR model 3 F 4 
mass of 4021 MeV, we predict a strong width for this 
state of just 8.3 MeV. This is partly a result of the cen- 
trifugal barrier (both open modes, DD and DD*, are G- 
waves), and is also because there is essentially no phase 
space for D*D*. If the GI masses are more accurate we 
would no longer expect a very narrow 3 F4 state, since 
the D*D* partial width increases rapidly above thresh- 
old. At the GI mass of 4095 MeV the 3 F 4 is predicted 
to have a partial width to D*D* of 38 MeV, and a total 
width of 60 MeV. 

The single well-established mechanism that might be 
exploited to produce a IF state is an El radiative transi- 
tion from the ■0(4159) to the 3 F2 X2- This transition how- 
ever has a rather weak, model-dependent partial width of 
20-50 keV (Table E), and the 3 F 2 is predicted to be the 
broadest IF state, with a total width of about 160 MeV. 

If an effective production mechanism for IF states is 
identified, the El decays of these states can be used 
to populate the narrow ID states; the El transitions 
3 F 4 -> 3 D 3 , 3 F 3 -> 3 D 2 and X F 3 -> X D 2 are predicted 
to have rather large radiative widths of ca. 300-400 keV 
(see Table ED. 



H. 2F and 1G states 

To complete our study we have evaluated the spectrum 
and decays of all states in the 2F and 1G multiplets (Ta- 
bles 1X1 VI and IXV|) , although the lack of a clear experi- 
mental route to these states suggests that it may be quite 
difficult to test these predictions. 

The multiplets are again predicted to be nearly degen- 
erate in both the NR and GI potential models, with some- 
what lower masses predicted by the NR model. The mean 
2F and 1G masses are approximately 4350 MeV and 
4225 MeV respectively in the NR model, and 4425 MeV 
and 4320 MeV in the GI model. 



Assuming NR model masses, these states are predicted 
to have total widths that range from a minimum of 
58 MeV for the 3 G 5 5 (4214) to a maximum of 180 MeV 
for the 2 3 F 2 % 2 (4351). The individual decay modes do 
show some interesting features; although the dominant 
modes are usually of S+S type, some S+P modes are 
important in 2F decays. For example, the largest decay 
modes of the 2 3 F 3 xa(4352) and 2 3 F 2 X2(4351) are both 
S+P, with branching fractions of B(x 3 — » DD 2 ) = 30% 
and B(% 2 — > DDi) = 59% respectively. 



VI. SUMMARY AND CONCLUSIONS 

We have computed the spectrum, all allowed El and 
some Ml electromagnetic partial widths, and all allowed 
open-charm strong decay amplitudes of the 40 charmo- 
nium states expected to ca. 4.4 GeV. These were evalu- 
ated using two potential models, the quark model formu- 
lation of electromagnetic transition amplitudes, and the 
3 Po strong decay model. 

The predictions of the spectrum should be useful in 
the identification of new states, and in tests of the 
Lorentz structure of confinement and the nature of spin- 
dependent forces. The transition amplitudes will be use- 
ful in searches for new higher-mass cc states, as they sug- 
gest which final states should be populated preferentially 
by the decays of a given initial cc state, as well as predict- 
ing the often characteristic strong decay amplitude ratios 
within a given final state. These predictions may also be 
useful in distinguishing between conventional charmonia 
and exotica such as charmed meson molecules (perhaps 
including the X(3872)) and cc hybrids. These results all 
implicitly test the accuracy of the pure-cc assumption 
made in these models, and can serve as benchmarks for 
calculations that relax this assumption. 

Future experimental measurements of strong partial 
widths and decay amplitude ratios can provide very im- 
portant tests of strong decay models, the 3 Po model in 
particular. We have discussed many specific examples of 
these tests in the previous sections; here we will summa- 
rize some of our most important observations. 

The -0(3770) is the lightest known cc state above open- 
charm threshold. We noted that the relatively large e + e~ 
width of the 0(3770) is difficult to understand if it is 
assumed to be a pure 3 Di cc state, but this problem is 
solved if the 0(3770) has a significant admixture of the 
2 3 Si cc basis state. This mixing also solves the 0(3770) 
total width problem, since the mixing angle required to 
fit the e + e~ width is consistent with the value required to 
fit the 3 Po model DD strong width to experiment. This 
mixing angle has been discussed in earlier references (see 
for example Refs. |75ll7qP . although the observation that 
the same mixing angle resolves both the e + e _ and total 
width discrepancies has not been noted previously. 

The relative branching fractions of the -0(4040) to DD, 
DD* and D*D*, although often cited as anomalous, are 
naturally explained in the 3 Pq model by a node in the DD 
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decay amplitude. A very important test of this strong de- 
cay model and the nature of the 0(4040) follows from a 
simple measurement of the ratios of the three indepen- 
dent D*D* decay amplitudes. The 3 Po model predicts 
that the 5 Fi amplitude is zero, and the ratio of the two 
P-wave amplitudes is 5 Pi/ 1 Pi = — 2v5, independent of 
the radial wavefunction, provided that the ^(4040) is a 
pure S-wave cc state. Similarly, the D*D* amplitude ra- 
tios in 0(4159) decays test the decay model and the as- 
sumed 2 3 D 1 assignment for this state. The 5 F! D*D* 
amplitude is predicted to be the largest for a pure 2 3 Di 
-0(4159), and the P-wave amplitude ratio is predicted to 
be ^i^Pi = -1/V5, again independent of the radial 
wavefunction. Finally, we noted that it may be possible 
to identify a higher-L cc state (the 3 F2 Xi member of the 
F-wave cc multiplet) in El radiative 0(4159) decays. 

The V(4415) has ten open-flavor modes, and their 
branching fractions have never been determined. We 
predict that the largest is the unusual S+P combina- 
tion DDi. It is especially notable that strong decays 
of the 0(4415) may provide a novel production mech- 
anism for the enigmatic D* (2317), through the above- 
threshold decay 0(4415 + e) -> D*D* (2317), as shown in 
Fig. One may also produce higher-mass charmonium 
states from the 0(4415) through electromagnetic transi- 
tions. For example, 0(4415) El decays are predicted to 
produce 3P states with branching fractions at the per mil 
level. These states could then be identified through their 
subsequent open-charm strong decays. 

As discussed in Section V.C, two new states near 
3940 MeV have recently been reported by Belle. 2P xj 
states are natural assignments, as they have theoretical 
quark model masses ranging from 3852 to 3972 MeV and 
widths of 50-100 MeV, similar to the reported values. 
A measurement of the relative branching fractions to DD, 



DD* and D*D* and comparisons to our 3 Po decay model 
predictions should allow the determination of the quan- 
tum numbers of these states, and will show whether they 
are indeed consistent with 2P cc assignments. 

A better determination of the properties of conven- 
tional charmonium states at higher masses is important 
not only because of the improved understanding of QCD 
that will follow (especially aspects of confinement and 
strong decays), but also because qualitatively different 
types of resonances are expected at these masses. These 
new states include charmed meson molecules and char- 
monium hybrids, and the identification of these novel 
excitations will obviously be easier if the conventional 
charmonium spectrum is well established. 

The recent results from B factories and new programs 
at BES, CLEO and GSI have led to a resurgence of in- 
terest in the physics of charmonium. We argue that a 
detailed experimental investigation of the spectrum of ex- 
cited charmonium states and their decay properties will 
considerably improve our understanding of the nonper- 
turbative aspects of QCD. 
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TABLE II: S — » P El radiative transitions in the NR and GI potential models. The masses are taken from Table [I] we use 
the experimental masses (rounded "input" column) if known, and for the x Pi h c we assume a mass of 3525 MeV, which is the 
c.o.g. of the 3 Pj Xi states. Otherwise, theoretical values are used. 



Multiplets 


Initial meson 


Final meson 


E-., 
NR 


(MeV) 
GI 


Tthy 

NR 


(keV) 

GI 


r expt (keV) 


2S — > IP 


ip (2 J Si) 


X2(1 3 P2) 


128. 


128. 


38. 


24. 


27. ± 4. 






Xi(l 3 Pi) 


171. 


171. 


54. 


29. 


27. ± 3. 






Xo(l 3 Po) 


261. 


261. 


63. 


26. 


27. ± 3. 






ftcCl'Pl) 


111. 


119. 


49. 


36. 




3S — > 2P 


^(3 3 Si) 


X2(2 3 P 2 ) 


67. 


119. 


14. 


48. 








Xi(2 3 Pi) 


113. 


145. 


39. 


43. 








Xo(2 3 Po) 


184. 


180. 


54. 


22. 






^So) 


ft c (2 1 Pi) 


108. 


108. 


105. 


64. 




3S — > IP 


V>(3 3 Si) 


X2(1 3 P2) 


455. 


508. 


0.70 


12.7 








Xi(l 3 Pi) 


494. 


547. 


0.53 


0.85 








Xo(l 3 Po) 


577. 


628. 


0.27 


0.63 








^(l'Pi) 


485. 


511. 


9.1 


28. 




4S — > 3P 


^(4 3 Si) 


X2(3 3 P 2 ) 


97. 


112. 


68. 


66. 








Xi(3 3 Pi) 


142. 


131. 


126. 


54. 








Xo(3 3 Po) 


208. 


155. 


0.003 


25. 






7?c(4 1 S ) 




104. 


106. 


159. 


101. 




4S — > 2P 


V'(4 3 Si) 


X2(2 3 P 2 ) 


421. 


446. 


0.62 


15. 








Xi(2 3 Pi) 


423. 


469. 


0.49 


0.92 








Xo(2 3 Po) 


527. 


502. 


0.24 


0.39 






7?c(4 1 S ) 


/i c (2 1 Pi) 


427. 


444. 


10.1 


31.3 




4S -> IP 


V^(4 3 Si) 


X2( 3 P 2 ) 


775. 


804. 


0.61 


5.2 








Xi( 3 Pi) 


811. 


841. 


0.41 


0.53 








Xo( 3 Po) 


887. 


915. 


0.18 


0.13 






?7c(4 1 S ) 


ZlcCPl) 


782. 


808. 


5.2 


9.6 
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TABLE III: IP and 2P El radiative transitions (format as in Table ITU. 



Multiplets 


Initial meson 


Final meson 


NR 


(MeV) 
GI 


NR 


(keV) 
GI 


r expt (keV) 


IP — > IS 


X2(1 3 P2) 




429. 


429. 


424. 


313. 


426. ± 51. 




Xi(l 3 Pi) 




390. 


389. 


314. 


239. 


291. ± 48. 




Xo(l 3 Po) 




303. 


303. 


152. 


114. 


119. ± 19. 






^(^So) 


504. 


496. 


498. 


352. 




2P — » 2S 


X2(2 3 P 2 ) 


^ (2 Si) 


276. 


282. 


304. 


207. 






Xi(2 3 Pi) 




232. 


258. 


183. 


183. 






Xo(2 3 P ) 




162. 


223. 


64. 


135. 






M2 x Pi) 


?7c(2 1 So) 


285. 


305. 


280. 


218. 




2P — > IS 


X2(2 3 P 2 ) 


J/^(rSx) 


779. 


784. 


81. 


53. 






Xi(2 3 Pi) 




741. 


763. 


71. 


14. 






Xo(2 3 Po) 




681. 


733. 


56. 


1.3 






^(2^1) 


^(l'So) 


839. 


856. 


140. 


85. 




2P -» ID 


X2(2 3 P 2 ) 


^ 3 (l 3 D 3 ) 


163. 


128. 


88. 


29. 








^ 2 (l 3 D 2 ) 


168. 


139. 


17. 


5.6 








iA(l 3 Di) 


197. 


204. 


1.9 


1.0 






Xi(2 3 Pi) 


^2(1 3 D 2 ) 


123. 


113. 


35. 


18. 








^(l 3 Di) 


152. 


179. 


22. 


21. 






Xo(2 3 Po) 


^(l 3 Di) 


81. 


143. 


13. 


51. 






/i c (2 1 Pi) 


r?2c(l 1 D 2 ) 


133. 


117. 


60. 


27. 
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TABLE IV: 3P El radiative transitions (format as in Table HB. 



Multiplets 


Initial meson 


Final meson 


NR 


(MeV) 

GI 


Ithy 

NR 


(keV) 

GI 


r cxpt (keV) 


3P -> 3S 


X2(3 3 P 2 ) 


tf(3 s Si) 


268. 


231. 


509. 


199. 






Xi(3 3 Pi) 




225. 


212. 


303. 


181. 






Xo(3 3 P ) 




159. 


188. 


109. 


145. 






h c (3 l Pi) 




229. 


246. 


276. 


208. 





3P -» 2S 


X2(3 3 P 2 ) 


^'(2 3 Si) 


585. 


602. 


55. 


30. 




Xi(3 3 Pi) 




545. 


585. 


45. 


8.9 




Xo(3 3 Po) 




484. 


563. 


32. 


0.045 




h c {Z x P x ) 


r?c(2 1 S ) 


593. 


627. 


75. 


43. 


on . i o 

or — > Id 




J/lp(L bl) 


1 fl/1 Q 
lU4o. 


1 OAQ 


o4. 


1 Q 




Xi(3 3 Pi) 




1013. 


1048. 


31. 


2.2 




Xo(3 3 Po) 




960. 


1029. 


27. 


1.5 






^(l'So) 


1103. 


1131. 


72. 


38. 


or — > ID 


X2(.i P2) 


y>3(2 D3) 


1 A T 


110 

118. 


1 1 O 


C 1 

•!.) 1 . 






tit fO^TS \ 
W2\£ r>i) 


100. 


Iz 1 . 


Q1 
OX. 


O Q 

y . y 






lp(2 Di ) 


IKK 

loo. 


1/11 

141. 


O 1 

2.x 


U-M 




Xl(3 fi) 




1 Iz. 


1 HQ 

lUo. 


Oo. 


00. 






^(2 3 Di) 


111. 


121. 


19. 


15. 




Xo(3 3 Po) 


V(2 3 Dr) 


43. 


97. 


4.4 


35. 






r?2c(2 1 D 2 ) 


119. 


109. 


99. 


48. 


3P -> ID 


X2(3 3 P 2 ) 


i/> 3 (l 3 D 3 ) 


481. 


461. 


0.049 


6.8 






^ 2 (1 3 D 2 ) 


486. 


470. 


0.0091 


0.13 






^(l 3 Di) 


512. 


530. 


0.00071 


0.001 




Xi(3 3 Pi) 


V 2 (1 3 D 2 ) 


445. 


452. 


0.035 


4.6 






^(l 3 Di) 


472. 


512. 


0.014 


0.39 




Xo(3 3 Po) 


V>(l 3 Di) 


410. 


490. 


0.037 


9.7 




h c (3 X Pi) 


^(l'Da) 


453. 


454. 


0.16 


5.7 
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TABLE V: ID and 2D El radiative transitions (format as in Table ITU . 



Multiplets 


Initial meson 


Final meson 


NR 


(MeV) 
GI 


Tthy 

NR 


(keV) 
GI 


Texpt (keV) 


ID -> IP 


^ 3 (1 3 D 3 ) 


x 2 (i 3 P 2 ) 


242. 


282. 


272. 


296. 






</> 2 (l 3 D 2 ) 


x 2 (i 3 P 2 ) 


236. 


272. 


64. 


66. 








Xi(l 3 Pi) 


278. 


314. 


307. 


268. 






V'(l 3 Di) 


X2(1 3 P 2 ) 


208. 


208. 


4.9 


3.3 


< 330 (90% c.l.) f76. 871 






Xi(l 3 Pi) 


250. 


251. 


125. 


77. 


280 ± 100 f76, 871 






Xo(l 3 Po) 


338. 


338. 


403. 


213. 


320 ± 120 [76, 87] 




ft c a(l 1 D 2 ) 


M! lp i) 


264. 


307. 


339. 


344. 




2D -> 2P 


</> 3 (2 3 D 3 ) 


X 2 (2 3 P 2 ) 


190. 


231. 


239. 


272. 






i/> 2 (2 3 D 2 ) 


X2(2 3 P 2 ) 


182. 


223. 


52. 


65. 








Xi(2 3 Pi) 


226. 


247. 


298. 


225. 






V'(2 3 Di) 


Xa(2 3 P 2 ) 


183. 


210. 


5.9 


6.3 








Xi(2 3 Pi) 


227. 


234. 


168. 


114. 








Xo(2 3 Po) 


296. 


269. 


483. 


191. 






77 c2 (2 1 D 2 ) 


/i c (2 1 Pi) 


218. 


244. 


336. 


296. 




2D -> IP 


^3(2 3 D 3 ) 


x 2 (i 3 P 2 ) 


566. 


609. 


29. 


16. 






</> 2 (2 3 D 2 ) 


x 2 (i 3 P 2 ) 


558. 


602. 


7.1 


0.62 








Xi(l 3 Pi) 


597. 


640. 


26. 


23. 






V>(2 3 Di) 


x 2 (i 3 p 2 ) 


559. 


590. 


0.79 


0.027 








Xi(l 3 Pi) 


598. 


628. 


14. 


3.4 








Xo(l 3 Po) 


677. 


707. 


27. 


35. 






h c2 {2 1 B 2 ) 




585. 


634. 


40. 


25. 




2D -> IF 


^3(2 3 D 3 ) 


X4(1 3 F 4 ) 


143. 


120. 


66. 


26. 








X3(1 3 F 3 ) 


136. 


114. 


4.8 


1.9 








X 2 (1 3 F 2 ) 


136. 


123. 


14. 


0.055 






V' 2 (2 3 D 2 ) 


X3(1 3 F3) 


127. 


110. 


44. 


19. 








X 2 (1 3 F 2 ) 


127. 


114. 


5.6 


2.4 






</>(2 3 Di) 


Xa(l 3 F 2 ) 


128. 


101. 


51. 


17. 






r? c2 (2 1 D 2 ) 




130. 


112. 


54. 


23. 
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TABLE VI: IF El radiative transitions (format as in Table [TTJ. 



Multiplets Initial meson Final meson E 7 (MeV) r t h y (keV) r exp t (keV) 







INK 




Mr? 
INK 


Lit 


X4(l F 4 ) 


^ 3 (1 3 D 3 ) 


209. 


239. 


332. 


334. 


X3(1 3 F 3 ) 


^3(1 3 D 3 ) 


217. 


240. 


41. 


38. 




^ 2 (1 3 D 2 ) 


222. 


251. 


354. 


325. 


X2(1 3 F 2 ) 


V>3(1 3 D 3 ) 


217. 


236. 


1.6 


1.4 




V>2(1 3 D 2 ) 


222. 


246. 


62. 


69. 




V>(l 3 Di) 


251. 


309. 


475. 


541. 




^(l'Da) 


221. 


249. 


387. 


321. 
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TABLE VII: 2F El radiative transitions (format as in Table H]}. 



Multiplets 


Initial meson 


Final 


meson 


NR 


(MeV) 
GI 


Tthy 

NR 


(keV) 
GI 


r expt (keV) 


2F -> 2D 


X4(2 3 F 4 ) 




;2 3 d 3 ) 


177. 


203. 


307. 


297. 






Xa(2 3 F 3 ) 


1p3 


;2 3 d 3 ) 


181. 


204. 


36. 


35. 








1p2 


;2 3 d 2 ) 


190. 


213. 


334. 


289. 






X2(2 3 F 2 ) 


1p3 


;2 3 d 3 ) 


180. 


200. 


1.4 


1.4 








1p2 


;2 3 d 2 ) 


189. 


209. 


58. 


50. 










;2 3 Di) 


188. 


222. 


306. 


295. 






^3(2^3) 


r\ci 




188. 


211. 


362. 


284. 




2F -> ID 


X4(2 3 F 4 ) 




:i 3 d 3 ) 


508. 


539. 


20. 


8.6 






X3(2 3 F 3 ) 




;i 3 d 3 ) 


512. 


539. 


2.3 


0.14 










;i 3 d 2 ) 


517. 


549. 


19. 


11. 






X2(2 3 F 2 ) 


i^s 


;i 3 d 3 ) 


511. 


536. 


0.090 


0.003 








1p2 


;i 3 d 2 ) 


516. 


545. 


3.2 


0.39 










;i 3 Di) 


542. 


604. 


20. 


18. 






h c3 (2 1 F 3 ) 


r/c2 


;i 1 D a ) 


516. 


548. 


22. 


9.9 




2F -> 1G 


X4(2 3 F 4 ) 


1p5 


1 3 G 5 ) 


132. 


112. 


54. 


22. 








i>4 


1 3 G 4 ) 


118. 


104. 


2.0 


0.84 








1p3 


1 3 G 3 ) 


110. 


101. 


0.025 


0.011 






X3(2 3 F 3 ) 


1p4 


1 3 G 4 ) 


122. 


105. 


43. 


18. 








l(>3 


1 3 G 3 ) 


113. 


102. 


2.3 


1.0 






X2(2 3 F 2 ) 


l(>3 


1 3 G 3 ) 


113. 


98. 


36. 


16. 






^3(2^3) 


r] C 4 


1 X G 4 ) 


123. 


106. 


47. 


20. 





TABLE VIII: 1G El radiative transitions (format as in Table Hill. 



Multiplets Initial meson Final meson E 7 (MeV) r t h y (keV) r exp t (keV) 







TVTT? 

INK 


VjxL 


Mr? 
INK 


\jL 


l*s( 1G 5) 


X4(l F 4 ) 


189. 


212. 


373. 


355. 


i/> 4 (l 3 G 4 ) 


X4(1 3 F 4 ) 


202. 


219. 


29. 


25. 




X3(1 3 F 3 ) 


194. 


217. 


382. 


349. 


i/> 3 (l 3 G 3 ) 


X4(1 3 F 4 ) 


210. 


222. 


0.66 


0.52 




Xa(l 3 F 3 ) 


203. 


220. 


37. 


31. 




X2(1 3 F 2 ) 


203. 


225. 


425. 


366. 


77c4(l 1 G 4 ) 




194. 


217. 


407. 


374. 
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TABLE IX: Ml radiative partial widths. The assumed masses are experimental where known, and otherwise are the theoretical 
predictions (see Table HJ. One exception is the ?7 c (3 1 So), for which we assume a mass of 4011 MeV (the mass of the known 
^(4040) minus the theoretical 3S splitting). We give results for the NR model (both with and without the recoil factor of 
jo(kr/2)) and the GI model, which includes the recoil factor. 



Initial Multiplet Initial meson 


Final meson 




(MeV) 




Tthy (keV) 




Texpt (keV) 






INK 




INK 


Mr? c a \ 


r 1 t 




is j/v>(i 3 Si) 




116. 


115. 


2.9 


2.9 


2.4 


1.1 ±0.3 


2S V'(2 3 Si) 


77c(2 1 So) 


48. 


48. 


0.21 


0.21 


0.17 








639. 


638. 


4.6 


9.7 


9.6 


0.8 ±0.2 


'7c(2 1 So) 


J/^(l 3 Si) 


501. 


501. 


7.9 


2.5 


5.6 




3S ^(3 3 Si) 




29. 


35. 


0.046 


0.046 


0.067 








382. 


436. 


0.61 


1.8 


2.6 








922. 


967. 


3.5 


8.7 


9.0 






V'(2 3 Si) 


312. 


361. 


1.3 


0.22 


0.84 






J/^(l 3 Si) 


810. 


853. 


6.3 


1.9 


6.9 




IV h c (2 Fi) 


/i 3r> S 
X2(l P2) 


OCA 

360. 


380. 


0.0 1 1 


0.0(0 


nil 
0.11 






Xi(l 3 Pi) 


400. 


420. 


0.058 


0.13 


0.36 






Xo(l 3 Po) 


485. 


504. 


0.033 


0.21 


1.5 




X2(2 3 P 2 ) 


Mi'Pi) 


430. 


435. 


0.067 


0.90 


1.3 




x'i(2 3 Pi) 


hcC^Pi) 


388. 


412. 


0.050 


0.51 


0.045 




Xo(2 3 Po) 


fccC^Pl) 


321. 


379. 


0.029 


0.19 


0.50 
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TABLE X: Open-flavor strong decays, 3S and 4S states. The i/)(4040) and ^(4415) masses (boldface) are taken from experiment; 
for the remaining (unknown) masses we assume the theoretical NR potential model values of Table H] 



Meson 


State 


Mode 


Tth„ \r e -r„t~\ (MeV) 


Amps. (GeV" 1/2 ) 


t/j(4040) 


3 3 Si 


DD 


0.1 


Pi = —0.0052 






DD* 


oo 


3 Pi — — n ftQ^d 






D*D* 


33 


: Pi — +0 0338 
5 Pi — —0 1510 
5 Fi = 






D S D S 


7.8 


J Pi = +0.0518 






total 


74 [52 ± 10] 




») c (4043) 


3^0 


DD* 


47 


3 P = -0.1139 






D*D* 


33 


3 Pn — —0 1489 






total 


80 




^(4415) 


4 d Si 


DD 


0.4 


Pi = +0.0066 






DD* 


2.3 


Pi = +0.0177 






D*D* 


16 


x Pi = -0.0109 
Pi = +0.0487 
5 Fi = 






DDi 


31 


3 Si = 

3 Di = +0.0933 






DDI 


1.0 


3 Si = +0.0168 
3 Di = 






T~"\T\* 


23 


JJl = — U.Uool 








U.U 

u 


bi = — o. ( ■ 1U 

3 t\ r\ 

Di = 






D S D S 


1.3 


Pi = —0.0135 






D S D S 


2.6 


Pi = +0.0212 






D S D S 


0. / 


Pi = +0.0027 
Pi = —0.0119 

5 Fi = 






total 


78 [43 ± 15] 




») c (4384) 


4 1 So 


DD* 


6.3 


3 P = +0.0299 






D*D* 


14 


3 Po = +0.0473 






DDS 


11 


= +0.0497 






DD| 


24 


5 D = -0.1014 






D S D* 


2.2 


3 Po = +0.0201 






D S *D* 


2.2 


3 Po = -0.0231 






D S D S * 


0.6 


1 S = -0.0136 






total 


61 
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TABLE XI: Open-flavor strong decays, 2P and 3P states. The initial meson masses are predictions of the NR potential model, 
Table □ 



Meson 


State 


Mode 


T thy (MeV) 


Amps. (GeV~ 1/2 ) 


y 2 (3972) 


2 d P 2 


DD 


42 


X D 2 = +0.0992 






DD* 


37 


3 D 2 = -0.1172 






D S D S 


0.7 


^2 = +0.0202 






total 


80 




Xi (3925) 


2 a Pi 


DD* 


165 


3 Si = +0.2883 










3 Di = -0.0525 


Xo(3852) 


2 J Po 


DD 


30 


L S = +0.1025 


M3934) 


2'Pi 


DD* 


87 


3 Si = -0.1847 










3 Di = -0.0851 




3 il p 2 


DD 


8.0 


x Do — —0 0330 






DD* 


2.4 


3 T>2 = +0.0191 






D*D* 


24 


5 S2 = +0.0592 










1 D2 = +0.0107 










5 D2 = -0.0282 










5 G2 = 






DDi 


1.1 


3 P2 = +0.0240 










3 F? = +0 0105 






DDi 


12 


3 p 2 = -0.0915 










3 F 2 = 






D S D S 


0.8 


J D2 = +0.0115 






D,D* 


11 


3 D2 = -0.0474 






D*D* 


7.2 


5 S2 = -0.0266 










: D2 = +0.0145 










5 D2 = -0.0384 










5 G 2 = 






total 


66 




Yi(4271) 


3 3 Pi 


DD* 


6.8 


3 Si = -0.0337 










3 Di = +0.0011 






D*D* 


19 


5 Di = -0.0632 






DD*. 


0.1 


J Pi — —0 0062 






D S D* 


9.7 


3 Si = -0.0287 










3 Dx = -0.0385 






D*D* 


2.7 


5 Di = -0.0356 






total 


39 




Yn (42021 


3 3 Po 


DD 


0.5 


i So = -0.0091 






D*D* 


43 


1 S = -0.0267 










5 D = -0.0997 






D S D S 


6.8 


1 S = -0.0374 






total 


51 




hn (4279) 




DD* 


3.0 


3 Si — +0 0216 










3-pj — in rxrxAR 

1 J 1 — tU.UU^O 






D*D* 


22 


3 Si = +0.0456 










3 Di = -0.0487 






DD*, 


28 


x Pi = -0.0943 






D S D S * 


15 


3 Si = +0.0222 










3 Di = -0.0539 






D S *D S * 


7.5 


3 Si = -0.0464 










3 Di = -0.0327 






total 


75 
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TABLE XII: Open-flavor strong decays, ID and 2D states. The i/>(3770) and ^(4159) masses (boldface) are taken from 
experiment; for the remaining (unknown) masses we assume the theoretical NR potential model values of Table H] 



Meson 


State 


Mode 


Tthy [Tempi] (MeV) 


Amps. (GeV" 1/2 ) 


ifa ("38061 


1 3 D 3 


DD 


0.5 


X F 3 = +0.0150 


i/>(3770) 


l^Di 


DD 


43 [23.6 ± 2.7] 


i Pi = +0.1668 


7/>3(4lD/ ) 


2 JJ3 


JJD 




r 3 = +U.U00I 






DD* 


50 


3 F 3 = -0.0997 






U D 


( 


P3 = —0.1249 
X F 3 = +0.0218 

5 tt 1 r\ n /i to 

13 = —0.0478 
5 H 3 = 






D S D S 


5.7 


13 = +0.0358 






u s u s 


i 


3t7 n nonce 
r 3 = — U.UzUo 






total 






V'2(4158) 


2 D2 


DD* 


34 


P2 = +0.0121 
3 F 2 = -0.0822 






D*D* 


32 


5 P 2 = -0.0660 
£2 = —0.0685 






* 

D S D S 


26 


P2 = +0.0983 
F2 = —0.0149 






total 


(10 




V>(4159) 


2 3 Di 


DD 


16 


'Pi = -0.0522 






DD* 


0.4 


Pi = +0.0085 






D*D* 


35 


'Pi = +0.0489 
Pi = —0.0219 
5 Fi = -0.0845 






D s Ds 


8.0 


Pi = —0.0427 






D S D S * 


14 


3 Pi = +0.0733 






total 


74 [78 ± 20] 




?7c2(4158) 


2'D 2 


DD* 


50 


3 P 2 = -0.0099 
3 F 2 = -0.1007 






D*D* 


43 


3 P 2 = -0.0933 
3 F 2 = -0.0593 






D S D S * 


18 


3 P 2 = -0.0802 
3 F 2 = -0.0182 






total 


111 
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TABLE XIII: Open-flavor strong decays, IF states. The initial meson masses are predictions of the NR potential model, Table[I] 



Meson 


State 


Mode 


r,/,,, fMeV") 


Amps. (GeV" 1/2 ) 


-\/ , f/1091 A 


i r 4 


nn 

uu 


D.o 


u4 — TU.uoiy 








1 1 


3p H r\ 090/1 

u4 - — — U.UZU4 






p* p* 
jj jj 


U.Uo 


5pi pi 0009 

JJ4 — — U.UUyZ 




















5ri 00 1 n — 5 
Lj4 = — Z.O • 1U 










5t n 

14 — U 








n 09 
u.uz 


lp , in 009Q 






^ „ i „ 7 

total 


o.o 






1 a T7 
1 £3 


PP* 

JJJJ 


OQ 

oo 


JJ3 = +U.loo4 










= — u.uiyv 






JJ JJ 


n o 


5 ta n n 1 of? 
JJ3 = — U.Ulob 










Lj3 = — z.4 • 1U 






total 


Q /| 






1 r 2 


nn 


Q8 

• 'O 


in,, in. 1 /1QO 






Pi Pi* 


( 


3n,, 1 n 1 000 






P * P* 
JJ JJ 


U. 1 


02 = u 










ip» 1 n noQO 

JJ2 — +U.UU0U 










5 D 2 = -0.0061 










5 G 2 = -3.1 ■ 10~ 4 






D S D S 


5.9 


^2 = +0.0464 








161 




/i c3 (4026) 


l A F a 


DD* 


61 


^Da = -0.1313 










3 G 3 = -0.0220 






D*D* 


0.1 


3 D 3 = -0.0129 



3 G 3 = -1.3 ■ 10" 4 



total 61 
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TABLE XIV: Open-flavor strong decays, 2F states. The initial meson masses are predictions of the NR potential model, Table[I] 



Meson 


State 


Mode 


r ihu (MeV) 

- 1 - Lily Y » y 


Amps. (GeV" 1/2 ) 


•w / A Q A Q\ 


z r 4 


JJJJ 


1 o 
Iz 


a /"i in none 






JJJJ 


O 1 

ol 


3/"i n ncrn 

w = — u.uo/y 






JJ JJ 


O 1 

zl 


5p\ n noot: 
JJ4 — — U.U00D 










1 in nonn 

= +u.uzuy 










L14 — — U.U410 










5t n 
I4 = U 






JJJJi 


n r, 
U.O 


3tti in nnoo 
r 4 = +U.UUoz 










3n 1 n nni c; 
H4 = +U.UU10 






L)\Ji 


•"> n 
Z.U 


3 17 nni to 
1* 4 = — U.U1 1 Z 










•3 T T fi 






JJJJ 2 


U.U4 


5t? n nntt 1 
r 4 — — U.UUol 










3tt 11 in - 4 
II4 = —1.1 ■ 1U 






"Pi* Tt* 


U.UO 


3 17 n nnc:Q 
1 4 = — U.UUoo 










3tt n 

XI4 — U 








o . u 


U4 — -hU.UZoO 






"n ~r\* 

u s u s 


4.0 


3ri n nonn 

Vj4 — — U.UZyU 






T\* Ft* 


i i 
1 1 


5t\ n ncco 

JJ4 — — U.U0O2 










In in nno/i 

Lj4 — +U.UUo4 










L14 — — U.UUOo 










5t n 

I4 — u 






ZOZCLl 


o i 




, (A QCO\ 


Z *3 


"Pi Pi* 

JJJJ 


07 

Z f 


JJ3 = — U.UzUo 










3n ri nc:n7 
= — U.UOy/ 






p\ * "Pi* 

JJ JJ 


Z f 


5r\ n no 1 
JJ3 = —U.UzlZ 










Vj3 = — U.Udod 






P»P\ * 


n c 


It? n n 1 1 n 

1*3 = — u.uiiy 






JJJJl 


9 /I 

o-4 


3T7 1 n nooT 








1 

l.o 


3tp n m A *7 
T3 = — U.U14^ 








Q O 
OZ 


5t> 1 n 1 /i 00 
J^3 = +U.14.5Z 










5t? 1 n nnoc: 
r 3 = +U.UUZ0 










5xj„ in 1 n— 4 

1 1 3 — — i.y ■ iu 






T>*r»* 

JJ JJ 


n 9 
u.z 


r3 ■ — — U.UUyi 






P» T\* 


lo 


JJ3 = +U.U4o4 










3/"i n no^n 

Ij3 — — U.UZOU 






JJ S JJ S 




1 J 3 — — U.U0Z0 










{j3 — — U.Ullz 






J^sJJsO 


U.U4 


It? n fiAQO 

r 3 — — U.UUoy 






tOXCLl 


1 1 n 

1 1U 




X2(4<351) 




JJJJ 


lo 


1 Ti n n^i /i 
JJ2 = — U.U440 






JJJJ 


o n 
Z.U 


3r~\ n ni ti 
JJ2 = — U.U1 / 1 






P\* Pt* 
JJ JJ 


1 1 

41 


n 

02 = u 










In 1 n ni 07 
JJ2 = +U.U1Z / 










5t~\ n nnnfi 

JJ2 = — u.uuyo 










5n n nooQ 
O2 — —U.Uozo 






P\P\ 
JJJJl 


lUo 


3r> 1 n onno 

j r , 2 = +u.zuyo 










3i? 1 n ni fiT 
r 2 — +U.U10 / 






TYlV 

JJJJ X 


U.o 


3r> n 










3i7 n ni 1 q 
r 2 = — U.Ullo 






P»P\ * 


r r 
0.0 


5t> 1 n ncno 

1^2 = +u.uoyz 






D*DS 


0.2 


3 F 2 = +0.0064 






D S D S 


1.1 


X D 2 = +0.0134 






D S D* 


6.9 


3 D 2 = +0.0365 






D*D* 


2.7 


5 S 2 = 










^2 = +0.0194 










5 D 2 = -0.0146 










5 G 2 = -0.0140 








180 
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TABLE XIV: Open-flavor strong decays, 2F states (cont.). 



Meson 


state 


Mode 


tthy (MeV) 


A , //it; — l/2\ 

Amps. (Gev ' J 


M4350) 


2 X F 3 


DD* 


34 


3 D 3 = +0.0173 
3 G 3 = -0.0689 






D*D* 


24 


3 D 3 = -0.0267 
3 G 3 = -0.0584 






DDS 


11 


^3 = -0.0516 






DDi 


0.02 


3 F 3 = -0.0026 






dd; 


0.03 


3 F 3 = -0.0035 






DD£ 


22 


5 P 3 = -0.1209 
5 F 3 = -0.0070 
5 H 3 = -1.8 ■ 10" 4 






D*D5 


0.006 


3 F 3 = -0.0018 






D S D S * 


12 


3 D 3 = -0.0379 
3 G 3 = -0.0297 






D*D* 


6.0 


3 D 3 = -0.0400 
3 G 3 = -0.0098 






DsD* 


0.3 


4 F 3 = -0.0108 






total 


109 
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TABLE XV: Open-flavor strong decays, 1G states. The initial meson masses are predictions of the NR potential model, TableQ 



Meson 


otate 


iviocie 


1 thy (Me V ) 


Amps. V ) 


^5 (4214) 


1 3 G 5 


DD 


10 


"Hs = +0.0397 






DD* 


6.4 


3 H 5 = -0.0342 






D*D* 


41 


5 F 5 = -0.0982 










1 H 5 = +0.0045 










5 H 5 = -0.0083 










5 Js = 






D S D S 


0.4 


^5 = +0.0086 






D S D S * 


0.0 


3 H 5 = -0.0027 






total 


58 




^4(4228) 


1 3 G 4 


DD* 


88 


3 F 4 = +0.1209 










3 H 4 = -0.0334 






D*D* 


19 


5 F 4 = -0.0643 










5 H 4 = -0.0153 






D S D S * 


3.5 


3 F 4 = +0.0307 










3 H 4 = -0.0031 






D S *D S * 


0.0 


5 F 4 = -1.5 ■ 10 -4 










5 H 4 = -6.0 • 10" 7 






total 


111 




V>a(4237) 


1 3 G 3 


DD 


63 


X F 3 = +0.0974 






DD* 


66 


3 F 3 = +0.1077 






D*D* 


13 


^3 = +0.0358 










5 F 3 = -0.0326 










5 H 3 = -0.0217 






D S D S 


10 


1 F 3 = +0.0443 






D S D S * 


3 


3 F 3 = +0.0289 






D S *D S * 


0.0 


!F 3 = +4.6 • 10~ 4 










5 F 3 = -4.2 • 10~ 4 










5 H 3 = -1.4 ■ 10~ 5 






total 


155 




^4 (4225) 


1 1 G 4 


DD* 


72 


3 F 4 = -0.1076 










3 H 4 = -0.0367 






D*D* 


24 


3 F 4 = -0.0732 










3 H 4 = -0.0136 






D S D S * 


3 


3 F 4 = -0.0268 










3 H 4 = -0.0033 






D S *D S * 


0.0 


3 F 4 = -2.2 • 10 -5 










3 H 4 = -1.7 • 10~ 8 






total 


99 





